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Abstract. Using the Riemann-Liouville and Caputo Fractional Standard Maps 
(FSM) and the Fractional Dissipative Standard Map (FDSM) as examples, 
we investigate types of solutions of non-linear fractional differential equations. 
They include periodic sinks, attracting slow diverging trajectories (ASDT), 
attracting accelerator mode trajectories (AMT), chaotic attractors, and cas- 
cade of bifurcations type trajectories (CBTT). New features discovered in- 
clude attractors which overlap, trajectories which intersect, and CBTTs. 
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1. Introduction 

In recent years fractional calculus (FC) and fractional differential equations (FDE) 
became very popular in many areas of science. The books dedicated to the appli- 
cations of FDEs published in 2010-2011 include [Tj[21|3] in physics in general, [H[S] 
in modeling and control, [B] in viscoelasticity, [7] in systems with long-range inter- 
action. A good review of applications of FC to chaos in Hamiltonian Systems is 
given in [8]. Because fractional derivatives are integro-differential operators, they 
are used to describe systems distributed in time and/or space: systems with long 
range interaction [TJ [HI HQ1 HH [T^l H3] , non-Mar kovian systems with memory ( |14) 
Ch.10, [HI [TU [T71 HH HI]), fractal media [2U], etc. Biological systems are prob- 
ably the best examples of systems with memory. As it has been shown recently 
[21\ 122] , even processing of external stimuli by individual neurons can be described 
by fractional differentiation. In some cases [23l [24j [25l [26] FDEs are equivalent to 
the Volterra integral equations of the second kind. This kind of equations (non 
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necessarily FDEs) is used in nonlinear viscoelasticity (see for example [2TJ [28] ) 
and in biology (for applications to mathematical models in population biology 
and epidemiology see [2"9l I3U] ). 

As in regular dynamics (in application to various areas in science), the non- 
linearity plays a significant role in fractional dynamics. Chaos and order in nonlin- 
ear systems with long range interactions were considered in J2HU1HUH2], nonlinear 
FDEs in application to the control were reviewed in books [4j [5] , nonlinear frac- 
tional reaction-diffusion systems considered in [3TJ [3H [33J [33] . Corresponding to 
the fact that in physical systems the transition from integer order time derivatives 
to fractional (of a lesser order) introduces additional damping similar in appear- 
ance to additional friction [T4J [35], phase spaces of the systems with fractional 
time derivatives demonstrate different kinds of structures similar to the attractors 
of the dissipative dynamical systems El [351 ISS] ■ 

As in the case of the systems which can be described by regular differential 
equations, general properties of systems described by FDEs can be demonstrated 
on the examples of maps which can be derived by integrating the FDEs (if it is 
possible) over a period of perturbation. Equations of some fractional maps (FM), 
which include among others fractional standard maps (FSM) and fractional dis- 
sipative standard maps (FDSM), were derived in a series of recent publications 
pj E3 US ES HE EH horn the corresponding FDEs. 

Two of the most studied examples in the regular case are the Chirikov stan- 
dard map (SM) EI] and the Zaslavsky dissipative standard map (DSM) [^2l B3] . 
The SM provides the simplest model of the universal generic area preserving 
map and the description of many physical systems and effects (Fermi accelera- 
tion, comet dynamics, etc.) can be reduced to the studying of the SM. The topics 
examined include fixed points, elementary structures of islands and a chaotic sea, 
and fractional kinetics E2E31E2]- Different properties of the DSM were discussed 
in [321 EH EE1 ESI (SHI and a rigorous proof of the existence of chaotic at- 
tractor in such type of systems was obtained in [5_H [33J. It was also proved in 
[52L 153] that the system of this type exhibits quasi-periodic attractors, periodic 
sinks, transient chaos with the dynamics attracting to the sink, and the existence 
of the SRB measure for the purely chaotic dynamics. All this stimulated the study 
of the FSMs [351 IS3] and FDSM [3j|] in hope to reveal the most general properties 
of the fractional dynamics. 

Two-dimensional maps investigated in [381 1391 154] were derived from the 
FDEs (see Section [2]) and, as a result, they are discrete maps with memory in 
which the present state of evolution depends on all past states. The earlier studies 
of maps with memory were done on one-dimensional maps which were not derived 
from differential equations [SSJ [351 IH1 (SHI (Ml HO] • Most results were obtained for 
the generalizations of the logistic map and the main general result is that the 
presence of memory makes systems more stable. The initial study of the FSM 
in [38j[39l[54] was concentrated on the investigation of the fixed/periodic points' 
stability of the FSM, new properties and new types of attractors. 
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The results obtained in [38l |39j [54j revealed new unusual properties of the 
fractional attractors of the FSM, which are different not only from the properties 
of the fixed/periodic points of the non-dissipative systems (the SM in our case), 
but also from the properties of attractors of the regular (not fractional) dissipative 
systems (e.g. the DSM). The most unusual observed property of the FSMs is the 
existence (and persistence) of the new type of attractors - cascade of bifurcations 
type trajectories (CBTT). A cascade of bifurcations, when with the change in the 
value of a parameter a system undergoes a sequence of period-doubling bifurca- 
tions, is a well known pathway of a transition form order to chaos. In the case 
of the CBTT the period doubling occurs without a change in a system parame- 
ter and is an internal property of a system. In Section [3] we summarize results of 
the extensive numerical investigation of three fractional maps to demonstrate the 
properties of fractional attractors. An analysis of the conditions for the CBTTs' 
appearance will be presented in the following publications. 

2. Basic Equations 

The standard map in the form 



p n+1 = p n - Ksinx n , x n+1 =x n +p n+1 (mod 2ir) (2.1) 
can be derived from the differential equation 

oo 

x + K S m(x)J2 6 {f~ in + £) ) = °' (2 ' 2) 

n=0 

where e — > 0+, following the steps proposed by Tarasov in 

The equations for the Riemann-Liouville FSM (FSMRL) were obtained in 
[3"T] and [55]. Following the steps proposed [2B], the FSMRL in the form which 
converges to the SM as a —> 2 can be derived from the differential equation with 
the Riemann-Liouville fractional derivative describing a kicked system 

oo 

Q D?x + Ksm{x)J2 S (f ~ ( n + £ )) = °> (l<a<2), (2.3) 

where e — ¥ 0+, with the initial conditions 

( £>rV)(0+)=Pi, { Dr 2 x)(Q+) = b, (2.4) 

where 

D?x(t) = D? ir a x(t) = 

— r— / -, V rr (n-Ka<n), (2.5) 

T(n-a)dt n J (t - t)"-^ 1 v ~ ; ' y ' 

Df = d n /dt n , and qI" is a fractional integral. 

After integration of equation (2.3) the FSMRL can be written in the form 

Pn+i =Pn - Ksinx n , (2.6) 
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t-ra+1 



r(a) ^ 



Pi+ 



(mod 2tt), (2.7) 



where 



V*{m) = m 
and momentum p(t) is denned as 

P(t) = 



i-k 



{m-l) a - k (2.8) 

oU? _1 a(t). (2-9) 

Here it is assumed that T = 1 and 1 < a < 2. The condition b = is required 
in order to have solutions bounded at t — for a < 2 [38 . In this form the 
FSMRL equations in the limiting case a = 2 coincide with the equations for the 
standard map under the condition xq — 0. For consistency and in order to compare 
corresponding results for all three maps (the SM, the FSMRL, and the Caputo 
FSM (FSMC)) all trajectories considered in this article have the initial condition 
x = 0. 

Following the steps proposed [35] , the FSMC in the form which converges to 
the SM as a — > 2 can be derived from the differential equation similar to (2.3| but 
with the Caputo fractional derivative 



'D?x + Ksm(x) 



0, 



(1< a< 2) 



where e — > 0- 



where 



with the initial conditions 
p(0) = OgDlxHO) = (Djx)(0) = po, x(0) = x , 



(2.10) 



(2.11) 



$D?x{t) = I 



n—a 
t 



D^x{r)dT 
T(n - a) J (t-r) Q -™ +1 



D?x(t) - 
(n — 1 < a < n) 



(2.12) 

Integrating equation (2.101 with the momentum defined as p = x and assum- 



ing T = 1 and 1 < a < 2, one can derive the FSMC in the form 

—V 

(a - 1) U 



Pn 



Pn 



r(c 



n-l 



1) sinx,- 



sina;„ , (mod 27r), 



a;«+i = x n + po 



r(a) ^ 

v ' i=0 



V^(n — i + 1) sinxj, (mod 27r). 



(2.13) 



(2.14) 



It is important to note that the FSMC ((2.13), ( |2.14 1) can be considered on a torus 
(x and p mod 2ir), a cylinder (x mod 2tt), or in an unbounded phase space, whereas 
the FSMRL (( |2.6[ ), (2.7)) can be considered only in a cylindrical or an unbounded 
phase space. The FSMRL has no periodicity in p and cannot be considered on a 



torus. This fact is related to the definition of momentum ( 2.9 1 and initial conditions 
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(2.4). The comparison of the phase portraits of two FSMs is still possible if we 
compare the values of the x coordinates on the trajectories corresponding to the 
same values of the maps' parameters. 
The DSM in the form 



Xn+l — X n + P, 



n+1, 



Pn+i = -bP n - Z sin(X„) 



(2.15) 



(2.16) 



can be derived integrating the differential equation of the the kicked damped ro- 
tator (see for example [55] ) 



X + qX = esm(X)^T5(t 



(2.17) 



Two forms of the FDSM were derived in [7J [371 132] ■ The form which has been 
investigated numerically in [39] 



_i » 

X n+ i = V P k +iW a (q, k-n-1), 

1 (a — 1) f— i 



k=0 



P n +i = -bP n -Zsm(X n ), 
where functions W a (a, b) are defined by 

W a {a, b) = a 1 - a e a{b+1) \v{a - 1, ab) - T(a - 1, a(b + 1)) 

and r(a, b) is the incomplete Gamma function 

/>00 

T(a,b) = / y^e-vdy 



has been derived from the following fractional generalization of (2.17) 

oo 

D?X(t) - q D?X(t) = esin(X) £ 6(t - n), 

where 



(2.18) 
(2.19) 

(2.20) 
(2.21) 

(2.22) 



n=0 



gel, 1 < a < 2, fi = a-l. 



Here 



H = (1 - e~ g )/q, Z = - l iee q . 
Further in this paper we will also use K — eexp(g) and V = —q. 



(2.23) 
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3. Fractional Attractors 

3.1. Standard Map: Fixed and Periodic Points 

It appears that the dependence of the SM's fixed point (0,0) stability properties 
on the map parameter K plays an important role when transition to the FSMs is 
considered. (0,0) SM fixed point is stable for < K < 4. At K = 4 it becomes 
unstable but period two (T = 2) antisymmetric trajectory 

Pn+l = -Pn, X n+ i = -X„ (3.1) 

appears which is stable for 4 < K < 2ir. At K = 2ir this T=2 trajectory becomes 
unstable but it gives birth to two T = 2 trajectories [4T1 14~1] with 

Pn+l = —Pn, x n+l = %n ~ T- (3-2) 

These T — 2 trajectories become unstable when K ps 6.59 (see Figure^), at 
the point where T — 4 stable trajectories are born. This period doubling cascade 
of bifurcations sequence continues with T = 8 and T = 16 stable trajectories 
appearing at K « 6.63 and K w 6.6344 correspondingly. It continues until at K w 
6.6345 all periodic points become unstable and corresponding islands disappear. 
This scenario of the elliptic-hyperbolic point transitions with the births of the 
double periodicity islands inside the original island has been investigated in (61) 
and applied to investigate the SM stochasticity at low values of the parameter 
K < 4. 

3.2. Phase Space at Low K (Stable (0,0) Fixed Point) 

Stability of the fixed point (0,0) has been considered in [311 [S3]. It has been shown 
that for fractional values of a (we consider only 1 < a < 2) this point turns from 
elliptic (for a = 2) into a sink which is stable for 

A'<K cl (a) = ^^, (3.3) 

where 

oo 

Vai = J2(-l) k+1 VZ(k) (3.4) 
fc=l 

and can be calculated numerically (see Figure|4^). It can be shown that AT cl (2) = 4, 
which corresponds to the SM case. The structure of the FSM phase space for K < 
K c i preserves some features which exist in the a = 2 case. Namely, stable higher 
period points, which exist in the SM case, still exist in the FSM, but they exist in 
the asymptotic sense and they transform into sinks or, in the FSMRL case, into 
attracting slow (p n ~ n 2 ~ a ) diverging trajectories (ASDT) (see Figure [lj. Pure 
chaotic trajectories disappear. Each attractor has its own basin of attraction. The 
traces of the SM chaotic sea exist in the following sense: initially close trajectories, 
which do not start from any basin of attraction, may fall into absolutely different 
attractors. In most of the cases attracting points themselves do not belong to 
their own basins of attraction: trajectory which starts from an attracting point 
may fall into a different attractor. The FSMRL trajectories which converge to 
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Figure 1. (a). 1000 iterations on each of 25 trajectories for the 
SM with K = 0.6. The main features are (0,0) fixed point and 
T=2 and T=3 trajectories; (b). 400 iterations on each trajectory 
with po = 2 + O.Q4i, < i < 50 for the FSMRL case K = 0.6, 
a = 1.9. Trajectories converging to the fixed point and ASDTs 
of period 2 and 3 are present; (c). 100 iterations on each FSMC 
trajectory with po — —3.14 + 0.0314i, < i < 200 for the same 
case as in Fig. [TJd (K — 0.6, a — 1.9) considered on a torus. In 
this case all trajectories converge to the fixed point, period two 
and period three stable attracting points. 



the fixed point follow two different routs. Trajectories which start from the basin 
of attraction converge according to the power law x n ~ n~ x ~ a and p n ~ n~ a , 
while those starting from the outside of the basin of attraction are attracting 
slow converging trajectories (ASCT) with x n ~ n~ a and p n ~ All FSMC 

trajectories converging to (0, 0) fixed point follow the slow power law x n ~ n x ~ a 
and p n ~ n 1 ^" (see Figure J2J). This convergence may correspond to the recently 
introduced for the fractional dynamic systems notion of the generalized Mittag- 
Leffier stability [62], for which the power-law stability is a special case. 

3.3. Phase Space at K c \ < K < K C 2 (Stable T = 2 Antisymmetric Trajectory) 

Numerical simulations confirm that, as in the SM case, in the FSMs a period T = 2 
antisymmetric trajectory (sink) exists (asymptotically) for K\(a) < K < Ki(a) 
(4 < K < 2ir in the SM case) with 

Pn+l = —Pn, x n+l = ~X n . (3.5) 

Asymptotic existence and stability of this sink is a result of the gradual transfor- 
mation of the SM's elliptic point with the decrease in the order of derivative from 
a = 2 (see Figure [3]). As in the fixed point case, there are two types of convergence 
of the trajectories to the FSMRL T — 2 sink: fast with 

Sx n ~ n~ l ~ a , Sp n ~ rT a (3.6) 

and slow with 

6x n ~n- a , Spn-^n 1 -" (3.7) 
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Figure 2. Convergence of trajectories to the fixed point (0,0) 
at K < Kef. (a). Two trajectories for the FSMRL with K = 2, 
a = 1.4, and 10 5 iterations on each trajectory. The bottom one 
withpo — 0.3 is a fast converging trajectory. The upper trajectory 
withpo = 5.3 is an example of the FSMRL 's ASCT; b). Evolution 
of the FSMC trajectories with p = 1.6 + 0.002i, < i < 50 for 
the case K = 3, a — 1.9. The line segments correspond to the nth 
iteration on the set of trajectories with close initial conditions; 
c). x and p time dependence for the FSMC with K = 2, a = 1.4, 
xq = 0, and po = 0.3. 




Figure 3. Stable antisymmetric x n +i = — x n , p n +i = —p n pe- 
riod T — 2 trajectories for K = 4.5: (a). 1000 iterations on each 
of 25 trajectories for the SM with K = 4.5. The only feature is a 
system of two islands associated with the period two elliptic point; 
(b). The FSMRL stable T = 2 antisymmetric sink for a = 1.8. 
500 iterations on each of 25 trajectories: p a — 0.0001 + 0.08i, 
< i < 25. Slow and fast converging trajectories; (c). The FSMC 
stable T — 2 antisymmetric sink for a = 1.8. 1000 iterations on 
each of 10 trajectories: p = -3.1415 + 0.628i, < i < 10. 

and 

Sxn-n 1 - 01 , Spn-n 1 - 01 (3.8) 
convergence in the FSMC case [S3]. 
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Let's consider the FSMRL and assume the asymptotic existence of the T = 2 
sink (xi,pi) and {—xi, —pi). In this case the limit n 

K 



Pi 



■ sm(x;), 



oo in (2.6) and (2.7) gives 
(3.9) 



K 



OO 



^ ' fe=i 



The equation for xi takes the form 



xi 



K 



-Vai sm(xi), 



2T(a) 

from which the condition of the existence of the T = 2 sink is 

2r(a) 



K > K cl (a) 



V al 



From (3.12) it follows that K^a) = K cl (a) 



(3.10) 



(3.11) 



(3.12) 



• Remark 1: The T = 2 sink exists only in the asymptotic sense. On a trajectory 
which starts from point (xo, po)=(xi, pi) X2 ^ x\ and P2 ^ Pi- 

• Remark 2: In spite of the fact that in the derivation of (3.9) and (3.10) we 
used the limit n — > oo, computer simulations demonstrate convergence of 
all trajectories to the limiting values in the very good agreement with (3.6), 



(3.7), and (3.1 



Remark 3: Direct computations show that Ki(a) 
FSMC too. 



K c i{a) is valid for the 



3.4. Phase Space at K > K c2 

The SM's T — 2 antisymmetric trajectory becomes unstable when K = 2n : at 
the point in phase space where a pair of T — 2 trajectories with x„+i = x n — n, 
Pn+i = —p n appears. Numerical simulations show (see Figure [5]) that the FSMs 
demonstrate similar behavior. Let's assume that the FSMRL equations (2.6) and 
(2.7) have an asymptotic solution 

p n = (-!)>, x n = xi-^-[l-(-l) n }. (3.13) 



Then it follows from (2.6) that relationship (3.9) pi = K/2siri(xi) holds in this 
case too. Simulations similar to those presented in Figure [2)3 |63| show that for 
K > K C 2 (see Figure |4k) the FSMRL has an asymptotic solution 



Pn = (-!)>,+ An 1 



(3.14) 



with the same A for both even and odd values of n. Substituting (3.14) in (2.7) 
and considering limit n — > 00 one can derive (see |63j ) 

7rr(o;) 



sin(x ; ) 



KVca 



(3.15) 
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Figure 4. 

obtained by the numerical simulations of equations ( |3.12| , ( |3.16[ ) , 
(3.11 1, (3.15), (|3.9|) and confirmed by the direct simulations of 



Critical values in the FSMs' (a, K)-space. The results 



the FSMs: (a). The (0,0) fixed point is stable below K = K cl 
curve. It becomes unstable at K — K c \ and gives birth to the 
antisymmetric T = 2 sink which is stable at K c \ < K < K c i. 
A pair of T = 2 sinks with x n +\ = x n — n, p n +i = ~Pn is 
stable in the band above K = K C 2 curve. Cascade of bifurcations 
type trajectories (CBTTs) are the only phase space features that 
appear and exist in the narrow band which ends at the cusp on 
the figure. (x c ,p c ) is the point at which the SM's T = 2 elliptic 
points with x n +± = x n — n, p n +i = ~Pn become unstable and 
bifurcate (see Sec. 3.1); (b). a-dependence of the T = 2 sinks' 
x-coordinate for K = 4.5. The transition from the antisymmetric 
sink (upper curve) to x n+ i = x n — n, p n +i — —p n point takes 
place at a = 1.73. Solid lines represent areas of stability; (c). 
K = 4.5 a-dependence of the T = 2 sinks' p-coordinatc. 



which has solutions for 

K > K c2 = ^ (3.16) 

Val 

(see Figure |4| and the value of A can also be calculated: 

These results are in a good agreement with the direct FSMRL numerical simula- 
tions and hold for the FSMC. 

3.5. Cascade of Bifurcations Type Trajectories 

In the SM further increase in K at K ss 6.59 causes an elliptic-hyperbolic point 
transition when T = 2 points become unstable and stable T = 4 elliptic points 
appear. A period doubling cascade of bifurcations leads to the disappearance of 
the corresponding islands of stability in the chaotic sea at K ps 6.6344. The cusp 
in Figure [4^, points to approximately this spot (a = 2, K w 6.63). Inside of the 
band leading to the cusp a new type of the FSM's attractors appears: cascade of 
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Figure 5. Stable x n+1 — x n — 7T, p n +i = — Pn period T = 2 
trajectories for X > if c2 : (a). 500 iterations on each of 50 tra- 
jectories for the SM with K = 6.4. The main features are two 
accelerator mode sticky islands around points (—1.379,0) and 
(1.379,0) which define the dynamics. Additional features - dark 
spots at the top and the bottom of the figure (which are clear 
on a zoom) - two systems of T — 2 tiny islands associated 
with two T — 2 elliptic points: (1.379, 7r), (1.379 — tt,— tt) and 
(tt - 1.379, tt), (-1.379, -tt); (b). Two FSMRL's stable T = 2 
sinks for K — 4.5, a — 1.71. 500 iterations on each of 25 trajec- 
tories: po = 0.0001 -I- 0.08i, < i < 25; (c). Two FSMC's stable 
T = 2 sinks for K = 4.5, a = 1.71. 1000 iterations on each of 10 
trajectories: p = -3.1415 + 0.628«, < i < 10. 



bifurcations type trajectories (CBTT) (see Figure [6|. In CBTTs period doubling 
cascade of bifurcations occurs on a single trajectory without any change in the map 
parameter. In some cases CBTTs behave like sticky islands of the Hamiltonian 
dynamics: occasionally a trajectory enters a CBTT and then leaves it entering 
the chaotic sea. Near the cusp CBTTs are barely distinguishable - the relative 
time trajectories spend in CBTTs is small. With decrease in a the relative time 
trajectories spend in CBTTs increases and at a close to one a trajectory enters a 
CBTT after a few iterations and stays there over the longest computational time 
we were able to run our codes - 500000 iterations. A typical FSMC's CBTT of this 
kind is presented in the Figure [Tj;. 

CBTTs, easily recognizable in the phase space of the FSMRL (Figure 6i), 
often almost impossible to distinguish in the phase space of the FSMC (Figure 6 d). 
But they clearly reveal themselves on a coordinate versus time (step of iteration 
n) plot (Figure [6J:). 

3.6. More FSM's attractors 

The main feature of the SM's phase space at K > 2ir, which defines dynamics and 
transport and is well investigated (see for example [64 ) , is the presence of the sticky 
accelerator mode islands. Slight decrease in the order of derivative turns accelerator 
mode islands into attracting accelerator mode trajectories (AMT) (Figure [7b 
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Figure 6. Cascade of bifurcations type trajectories (CBTT) 
from iterations on a single trajectory with K — 4.5, a — 1.65, xq — 
0, and po = 0.3: (a). 25000 iterations for the FSMRL case. The 
trajectory occasionally sticks to one of the cascade of bifurcations 
type trajectories but later always returns to the chaotic sea; (b). 
30000 iterations for the FSMC case. The CBTTs can hardly be 
recognized on the full phase portrait but can be seen on the x of n 
dependence in Figure |6j:; (c). Time dependence of the coordinate 
x in Figure [6)3. 

Depending on the values of K and a trajectories either permanantly sink into this 
attractor, or intermittently stick to it and return to the chaotic sea. AMTs do not 
exist (we were unable to find them) for a < 1.995. In the (a, if)-plot (Figure [4^) 
the area of the AMT existence is a part (K > 2tt) of a very narrow strip near a = 2 
line. In the area of the a, if-plot above the CBTT strip (we considered only K < 7) 
different kinds of chaotic attractors (or possibly pure chaotic behavior) can be 
found. Two examples are presented in Figure[7j Figure]?}; contains two overlapping 
attractors one of which is a CBTT. Existence of the overlapping attractors and 
intersecting trajectories is a feature of maps with memory which is impossible in a 
regular dynamical systems. The simplest example of intersecting trajectories can 
be constructed by starting one trajectory at an arbitrary point (xq,Po) and then 
the second at (x' ,p' ) = (xi,pi). 

3.7. The Fractional Dissipative Standard Map (FDSM) 

Consideration of the FDSM in [39 demonstrated that already a small deviation 
of a from 2 leads to significant changes in properties of the map. For example, a 
window of ballistic motion which exists in Zaslavsky Map near K = 4tt is closing 
when a < 1.996. Further decrease in a at K m Air produces different kinds of 
chaotic attractors and sinks [39 . 

Some of the new features we observed, which are presented in Figure [HJ 
include inverse cascade of bifurcations and trajectories which start as a CBTT 
and then converge as an inverse cascade of bifurcations. In phase space those 
trajectories may look like chaotic attractors. 
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FIGURE 7. (a). A single FSMRL ballistic trajectory for K = 
1.999, a — 6.59; (b). Seven disjoint chaotic attractors for the case 
FSMRL with K = 4.5 and a = 1.02. 1000 iterations on each of 20 
trajectories: x = andp = 0.0001 + 1.65i, < % < 20; (c). 20000 
iterations on each of two overlapping independent attractors for 
the FSMC case with K = 4.5 and a = 1.02. The CBTT has 
Po = —1.8855 and the chaotic attractor po = —2.5135 (xq = 0). 




Figure 8. The FDSM (in all cases r = 5 and Q = 0): (a). 
Phase space for the case K = 9.1, a = 1.2; (b). 70000 iterations 
on a single trajectory for the same case as in Figure [8^, (K = 9.1, 
a = 1.2). A cascade of bifurcations appears from chaos and then 
converges into a periodic trajectory; (c). 70000 iterations on a 
single trajectory for the case K = 10, a = 1.2. A chaotic attractor 
appears in the FDSM's phase space after a chaotic trajectory 
through an inverse cascade of bifurcations forms a period T = 3 
trajectory. 

4. Conclusion 

We studied three different fractional maps which describe nonlinear systems with 
memory under periodic perturbations (kicks in our case). The types of solutions 
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which we found in all maps include periodic sinks, attracting slow diverging tra- 
jectories, attracting accelerator mode trajectories, chaotic attractors, and cascade 
of bifurcations type trajectories. New features discovered include attractors which 
overlap, trajectories which intersect, and CBTTs. 

The models with a similar behavior may include periodically kicked media 
which can be described by FDEs: viscoelastic materials [B] or dielectrics [55]. Ex- 
periments can be proposed to observe different kinds of solution, for example 
CBTTs, in those media. 

Nonlinear models with memory which have chaotic and periodic solutions, 
where cascades of bifurcations considered as a result of a change in a system 
parameter, are used in population biology and epidemiology [55]. Our calculations 
show that such behavior can be a consequence of the essential internal properties 
of the systems with memory and suggest construction of new fractional models of 
biological systems. 

Our computational results are based on the runs over more that 10 5 periods of 
perturbations. Due to the integro-differential nature of the fractional derivatives, 
at present time it is impossible to solve FDEs over more than 10 3 periods of 
oscillations. Our results suggest a new way of looking for a particular solutions of 
the FDEs. For example, we may suggest that CBTTs can be found at the values 
of system parameters, where in corresponding integer system a period doubling 
cascade of bifurcations leads to the disappearance of a system of islands. 
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